Thermodynamic properties of separable square-wave potentials 
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Exact analytic solutions to the Schrodinger equation for an electron moving in three dimensional 
potentials have been studied. These solutions can correspond to metals, semiconductors, or insu- 
lators. We show that there is an efficient method to calculate the electron density of states for 
this class of potentials. From the density of states, the temperature dependence of thermodynamic 
properties such as the chemical potential and the specific heat were determined. Ten thousand cubic 
separable potentials were considered. This data makes it possible to identify trends in how the form 
of the potential is related to the thermodynamic properties of a material. 

PACS numbers: 71.20.-b, 65.40.-b 
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I. INTRODUCTION 

To calculate the band structure of a crystal, some 
initial guess of the electron wave-functions is required. 
One of the rare cases where analytic solutions to the 
Schrodinger equation are known for a semiconductor or 
an insulator are separable potentials that can be written 
as the sum of three one-dimensional square-wave poten- 
tials. Such separable potentials were first considered by 
Kronig and Penney [l|. There is a simple analytic expres- 
sion for the energy - wave-number dispersion relation 
for the well-known problem of an electron moving in a 
one-dimensional square-wave potential. Here we show 
that there is also an analytic expression for the density 
of states for the one-dimensional problem. These results 
can be simply combined to produce the dispersion rela- 
tions and density of states of separable three dimensional 
potentials. A limiting case of this class of potentials is the 
widely used free electron model which occurs when the 
amplitude of the square-wave is zero. As the amplitude 
of the potential increases, the potentials can correspond 
to metals, semiconductors, or insulators. Normally de- 
termining the electronic density of states from the band 
structure is a computationally intensive process that re- 
quires sampling many points of the Brillouin zone. For 
the separable square- wave potentials, the electronic den- 
sity of states is however easily calculable as a convolu- 
tion of three analytically known functions. This makes 
it easy to calculate the thermodynamic properties corre- 
sponding to this class of potentials. From the density of 
states we can calculate the temperature dependence of 
thermodynamic quantities such as the chemical poten- 
tial, the internal energy, the specific heat, the entropy, 
and the Helmholtz free energy. Half a million plots of 
the thermodynamic properties of separable potentials are 
available as online supplementary material^. The col- 
lective data set allows us to consider how an increase in 
the amplitude of the potential causes a transition from a 
metal to an insulator. There are cases where increasing 



the amplitude of the potential causes a transition from 
a metal to an insulator, then back to a metal and then 
finally back to an insulator. 

While band gaps appear in the one-dimensional prob- 
lem for an arbitrarily small amplitude of the periodic 
potential, a finite amplitude is needed for the creation of 
band gaps in three dimensions. In section IHII we quan- 
tify how large the amplitude of the potential must be 
for bands to be formed for the specific class of three- 
dimensional potentials we have considered. This mini- 
mum amplitude is. 



V = 0.9- 



(1) 



Here V is the amplitude of the potential, h is the reduced 
Planck's constant, m is the mass of an electron, and a 
is the lattice constant. The amplitude of the potential 
necessary for the formation of bands depends strongly 
on the lattice constant. 



II. SOLUTIONS OF THE SCHRODINGER 
EQUATION FOR THE KRONIG-PENNEY 
POTENTIAL 

This section reviews the solutions to the Kronig- 
Penney model[^,|3] and derives an expression for the one- 
dimensional density of states which will be used in section 
mil to construct the dispersion relations and densities of 
states for three-dimensional separable potentials. Fig- 
ure [1] shows the one-dimensional potential that was first 
considered by Kronig and Penney [Ij. 

Because of the translational symmetry of the potential, 
the eigenfunctions of the Hamiltonian are simultaneously 
eigenfunctions of the translation operator. The eigen- 
functions of the translation operator T can be readily 
constructed from any two linearly independent solutions 
of the one-dimensional Schrodinger equation. A conve- 
nient choice is, 
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FIG. 1: One-dimensional square-wave potential. 



The solutions in region 1 (0 < x < 6) are, 

I I \ n \ I I \ sin(fcia;) 
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while the solutions in region 2 (& < a; < a) areQ, 

■01(2^) = cos(fc2(a; — 6)) cos(fci6) 
ki 

— —sm(k2(x — 6)) sm(kib), 
k2 
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where S = \J o? — \ and 



a = ^1 (a) H — ^ = 2 cos {k^ (a — 6) ) cos (fci 6) 
ax 

^-f ^)sin(fc2(a-6))sin(fci6). (8) 

fcl fc2/ 

If periodic boundary conditions are used for a potential 
with N unit cells, then applying the translation operator 
N times brings the function back to its original position 



The eigenvalues of the translation operator are therefore 
the solutions to the equation — 1. These solutions 
are. 



Aj = exp 



z27rj 



exp 
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j — exp {ikjo) , (10) 



■02(2^) = T~ cos(A;2(x — h)) sin(A:i6) 
ki 

1 



- sin(fc2(x - 6)) cos(fci6). (4b) 



Here k^ = ^J2m {E-Vr) /h?; r = 1,2 is the label of 
the region, and E is the energy. For energies where kr 
is imaginary, the solutions are still real since cos(i6') — 
cosh(0) and sm{i9) = isinh(6'). 

Any other solution can be written as a linear combi- 
nation of ipiix) and ip2ix). In particular, 0i(a; + a) and 
"02(2; + a) can be written in terms of ipii^) ^nd "02 (a^)- 
These solutions are related to each other by the matrix 
representation of the translation operator [6j. 



(5) 



The elements of the translation matrix can be determined 
by evaluating Eq. ([S]) and its derivative at x = 0. 
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where j is an integer between ~N/2 and N/2, L = Na 
is the length of the crystal, and kj = 2Trj/L are the al- 
lowed k values in the first Brillouin zone. The dispersion 
relation can be determined by first calculating a for a 
specific energy and then solving Eqs. ([7]) and for 
the wave- number [21, 



fc = ± — tan 

a 



V4-a2 



(11) 



The dispersion relation can be used to determine the 
density of states which is needed to calculate the ther- 
modynamic properties of a system of noninteracting 
electrons @. The one-dimensional density of states in 
k space is D{k) ^ 2/tt and thus the density of states in 
energy is, 
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The eigenfunctions and eigenvalues A of the translation 
operator are. 



where 
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dx 
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X± = -{a±S), (7) 
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FIG. 2: (a) The energy - wave-number dispersion relation. The dashed line is the Fermi energy, (b) The density of states, (c) 
The internal energy density (solid line) and Helmholtz free energy density (dashed line), (d) The chemical potential (solid line) 
and the specific heat (dashed line). All of the plots were drawn for a square- wave potential with the parameters: Vi = OeV, 
V2 = 12.5 eV, a = 2 X 10^^" m, 6 = 5 x 10^^^ m, and an electron density of n = 3 electrons/primitive cell. 
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Thermodynamic properties such as the chemical po- 
tential, the internal energy, the specific heat, the entropy, 
or the Helmholtz free energy of a system of noninteract- 
ing fermions can be calculated from the electron density 
n, the density of states D{E), and the temperature T. 

The chemical potential fi is implicitly defined by the 
relation, 



D{E)F{E)dE, 



where F{E) is the Fermi function, 

1 



(14) 
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Here fc^ is Boltzmann's constant. Once the chemical po- 
tential has been determined, it can be used to calculate 
the internal energy density u and the Helmholtz free en- 
ergy density /. 



= / ED{E)F{E)dE 



(16) 



f = Hn- ksT / D{E) In 1 + exp 
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(17) 

Finally, the specific heat and the entropy are given by 
the partial derivatives. 
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df 
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The band structure of a one dimensional potential and 
the corresponding thermodynamic properties are plotted 
in Fig. H 



III. THREE-DIMENSIONAL SEPARABLE 
POTENTIALS 

In this section, 5 of the calculated band structures 
that are available in the supplementary material are pre- 
sented. This illustrates the variety of results that can 
be obtained with this simple model. The five examples 
are a free electron g metal where the charge carriers 
at the Fermi surface are electron-like, a metal where the 
charge carriers at the Fermi surface are hole-like, a di- 
rect band gap semiconductor, and an indirect band gap 
semiconductor . 
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The dispersion relation and density of states for any 
three-dimensional potential of the form 



U3dix,y,z) = U^{x) + Uy{y) + U,{z) 



(19) 



are easily calculated from the one-dimensional results [9|. 
The energy of an electron in a three-dimensional sep- 
arable potential is the sum of the energies of the 
constituent one-dimensional potentials. The three- 
dimensional density of states is the convolution of the 
three one-dimensional densities of states. Once the three- 
dimensional density of states is known, the thermody- 
namic quantities can be calculated as outlined above. 

The one-dimensional bands can be indexed by integers, 
1 corresponding to the band with the lowest energy, 2 to 
the band with the next lowest energy, etc. The three- 
dimensional bands can then be indexed by the three in- 
tegers that correspond to the one-dimensional bands that 
make up the three-dimensional band. The three dimen- 
sional band with the lowest energy is the 111 band. For 
cubic crystals, the next three bands 211, 121, and 112 
are degenerate in energy. 

Figure [3] shows the band structure, the density of 
states, and corresponding thermodynamic properties for 
a constant potential that corresponds to a free electron 
gas. The band structure in Fig. [3^ is plotted along a 
path in the Brillouin zone going from the M point (0.5, 
0.5, 0) through the points F (0, 0, 0), X(0.5, 0, 0), M, 
R (0.5, 0.5, 0.5), X, F, ending at the point R. All of the 
standard results for a free electron gas are reproduced 
by the numerical calculation. The dispersion relation is 
parabolic; the density of states increases with the square 
root of the energy and the specific heat is a linear func- 
tion of the temperature. 

Figures 2] and [S] show the results for a potential with 
an amplitude of 4.2 eV, a = 0.2 nm, and b = 0.09 nm. 
Gaps open in the dispersion relation at the Brillouin 
zone boundaries and some kinks appear in the density 
of states. These kinks make both the analytical and nu- 
merical evaluation of the thermodynamic properties dif- 
ficult. Both Figs, m and [S] correspond to metals but in 
Fig. U the electron density is 2 electrons/primitive cell 
and the density of states is increasing at the Fermi en- 
ergy while in Fig. [5]the electron density is 1 electron/unit 
cell and the density of states is decreasing at the Fermi 
energy. The material in Fig. 2] has a Fermi surface that 
consists of electron-like states and the chemical potential 
decreases with increasing temperature as it does for the 
free electron gas of Fig. [3] The material in Fig. [S]has a 
Fermi surface that consists of hole-like states and conse- 
quently the chemical potential increases with increasing 
temperature. 

Figure [S] shows the band structure for a material with 
two band gaps. The lower band gap is indirect with a 
band gap energy of 3.9 eV. The upper band is direct 
with a band gap energy of 0.4 eV. If the material shown 
in Fig. inihad 2 electrons/primitive cell, the lowest band 
would be completely filled and it would be an insulator. 
For an electron density of 8 electrons/primitive cell, the 



Fermi energy lies in the second band gap (illustrated by 
the dashed line in Fig. [5^. The thermodynamic prop- 
erties in Figs. and [BJl were calculated assuming an 
electron density of 8 electrons/primitive cell. Since this 
material is a semiconductor, the electronic contribution 
to the specific heat is negligible at low temperature and 
increases exponentially at high temperatures. The chem- 
ical potential of a semiconductor is a linear function of 
the temperature. 

Figure [7] shows the band structure for a potential where 
the 311 band has moved down lower than the 221 band. 
In this case the lower band is indirect with a band gap 
of 4.7 eV and the upper band gap is also indirect with 
a band gap energy of 0.5 eV. For an electron density 
of 8 electrons/primitive cell, the Fermi energy lies in the 
second band gap where the dashed line is drawn and the 
material is an indirect band gap semiconductor. 

The dispersion relations, densities of states, and ther- 
modynamic properties like those shown in Fig. [3][7]were 
calculated for 10000 cubic potentials. For these calcula- 
tions, it is convenient to normalize the Schrodinger equa- 
tion so that length is measured in terms of the lattice 
constant a. When this is done, energies are measured 
in units of n^h? /2ma^ and there are only two indepen- 
dent parameters in the problem, which can be taken to 
be 6 = b/a and V = 2ma^V/n^h^. Here V = V2 - Vi. 
The calculations were performed for normalized param- 
eters in the range b — (0, 1) and V = (0, 27). It becomes 
more difficult to perform the numerical calculations of the 
thermodynamic properties for large values of V. The up- 
per limit of V — 27 corresponds to the amplitude where 
our calculations of the thermodynamic properties become 
unreliable. 

Figure [5] shows the sizes of the first two band gaps 
that occur. The three-dimensional potential must have a 
finite amplitude for there to be a band gap. 



IV. THERMODYNAMIC PROPERTIES OF THE 
METALS 

The band structures of the potentials that were studied 
always correspond to metals when there is no band gap. 
This region includes the free electron model a.t V = 
as well as (5-function potentials in the regions near b = 
and b — I. When band gaps occur, the situation is more 
complicated. The potentials can correspond to metals, 
semiconductors, or insulators depending on the electron 
density. We compared the thermodynamic properties of 
these metals to results obtained by using the Sommer- 
feld expansion. Sommerfeld showed that the electronic 
contribution to the thermodynamic properties of metals 
can be approximated in terms of just two quantities: the 
density of states at the Fermi energy D{Ep) and the 
derivative of the density of states at the Fermi energy [loj 



D'iEp) = 



dD{E) 



dE 



(20) 
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FIG. 3: The band structure and thermodynamic properties of a free electron gas. (a) The energy-momentum dispersion 
relation. The dashed line is the Fermi energy, (b) The density of states, (c) The internal energy density (solid line) and the 
Helmholtz free energy density (dashed line), (d) The chemical potential (solid line) and the specific heat (dashed line). All 
plots are for a cubic potential with the parameters: Vi = OeV, V2 ~ OeV, a — 0.2 nm, b = 0.09 nm, and an electron density of 
n = 2 electrons/primitive cell. 
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FIG. 4: The band structure and thermodynamic properties for a cubic potential Vi = eV, V2 ~ 4.2 eV, a — 0.2 nm, b — 0.09 nm 
and n — 2 electrons/primitive cell, (a)-(d) as in Fig. [S] 
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FIG. 5: The band structure and thermodynamic properties for a cubic potential Vi = eV, V2 ~ 4.2 eV, a = 0.2 nm, b — 0.09 nm 
and n — 1 electrons/primitive cell, (a)-(d) as in Fig. [S] 
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FIG. 6: The band structure and thermodynamic properties for a cubic potential Vi = OeV, V2 ~ 8.8 eV, a — 0.4 nm, b — 0.21 nm 
and n — 8 electrons/primitive cell. This band structure corresponds to a direct band gap semiconductor, (a)-(d) as in Fig. [S] 
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FIG. 7: The band structure and thermodynamic properties for a cubic potential Vi = eV, V2 = 9.5 eV, a = 0.4 nm, 6 = 0.19 nm 
and n = 8 electrons/primitive cell. This band structure corresponds to an indirect band gap semiconductor, (a)-(d) as in Fig. 
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FIG. 8: (a) The band gap of the lowest gap between band 111 and the degenerate bands 112, 121, and 211 in units of n^fi^ /2ma? . 
This band gap is always indirect, (b) The second band gap. The dashed lines in (b) indicate the transitions from direct to 
indirect band gaps. The indirect band gaps are in between the two dashed lines. Note that there is a reentrant regime at 
h — 0.3 where an increase in the amplitude V causes a transition from metal to insulator, then back to a metal and finally to 
an insulator. 



To lowest order in the temperature T, the Sommerfeld 
expression for the chemical potential, the internal energy, 
specific heat, entropy, and Helmholtz free energy are. 



D'{Ef) 



l{T = 0) 



TT^DjEF) 

6 



{kBTY 



(21) 



(22) 



(23) 



Van Hove singularities flT']. If the Fermi energy is within 
an energy fc^T of one of the kinks, then the assump- 
tions of the Sommerfeld theory are invalid. The numeri- 
cal simulations in the supplementary material show that 
the Sommerfeld approximation is valid in the majority 
of cases where there are no band gaps. When gaps are 
present, the bands are narrow and it is more likely that 
the Sommerfeld theory fails. In the case of extremely 
narrow bands it becomes difficult to determine the tem- 
perature dependence of the thermodynamic properties 
reliably because of the discontinuities in the density of 
states. 



(24) 



THERMODYNAMIC PROPERTIES OF THE 
SEMICONDUCTORS AND INSULATORS 



f~u(T^ 0) 



tt^D{Ef) 
6 



(kBTf 



(25) 



Figure [9] shows the value of the density of states at the 
Fermi energy and its derivative for an electron density 
of n = 1 electron/primitive cell. This figure makes it 
possible to estimate the thermodynamic properties of the 
metals by substituting the values for D{Ef) and D'{Ef) 
into Eqs. (EB)-®. 

The Sommerfeld expansion assumes that the density of 
states is a smoothly varying function in an energy range 
about kBT wide in the vicinity of the Fermi energy. How- 
ever, the density of states must contain kinks known as 



If the Fermi energy falls in a band gap so that the 
material is a semiconductor or an insulator, the density 
of states at the Fermi energy and its derivative are not 
suitable for describing the thermodynamic quantities. In- 
stead, it is common to specify the band gap and the ef- 
fective masses of the electrons and holes. 

For cubic crystals, the first gap appears above the 
111 band and beneath the next three bands (211, 121, 
and 112) which are degenerate due to symmetry reasons. 
Above these three degenerate bands, a second energy gap 
sometimes appears. The Fermi energy lies in the first gap 
for an electron density of 2 electrons per unit cell and it 
can only be in the second gap for an electron density of 
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FIG. 9: The density of states at the Fermi energy and the derivative of the density of states at the Fermi energy for an electron 
density of n = 1 electrons/primitive cell. It was not possible to determine the density of states at the Fermi energy reliably in 
the crosshatched regions. Similar plots for electron densities 2-10 electrons/primitive cell are available in the supplementary 
information online. 



8 electrons per unit cell. 

The first gap is always indirect with the maximum of 
the valence band occurring at R. The minimum in the 
conduction band for the first gap occurs at X. The second 
band is sometimes direct and sometimes indirect. The 
valence band maximum of the second gap is always at M. 
When the second band gap is direct, the conduction band 
minimum is also at M (see Fig. [6]) but when the second 
band gap is indirect, the conduction band minimum is at 
r (see Fig. [21). 

The effective masses near the top of the valence band 
and bottom of the conduction band can be found by lin- 
earizing a given by Eq. ([5]) near the band edges and 
inserting this into Eq. (fTTj) . The effective masses of elec- 
trons and holes are, 




and 




E=E^. 



E=E„ 



(26a) 



(26b) 



Here Ec is the energy at the bottom of the conduction 
band and Ey is the energy at the top of the valence band. 
The effective masses are plotted in Fig. [TOl 

The thermodynamic properties of semiconductors or 
insulators are typically calculated using the Boltzmann 
approximation. In this approximation, the density of 



states near the Fermi energy is described by the func- 
tion. 



'{2mlf'^^^ E<Ey 
D{E) = {0 E„<E<Er 



(27) 



(2m:) 



*^l3/2 ^E-Eg 



Ec<E 



In the Boltzmann approximation it is further assumed 
that in the conduction band the Fermi function in Eqs. 
([ni) . ([TB)) -([T^ can be replaced by a Boltzmann factor 
F{E) Ri exp{{fi — E)/kBT) while in the valence band 
F{E) « 1 - exp{{E-fj,) /ksT). The result for the chemi- 
cal potential of a semiconductor in this approximation is 
found in many textbooks [l^, [13 ■ The chemical potential 
of a semiconductor or insulator is a linear function of the 
temperature. 



Ey + Ec 



-|fc,rin(4 

4 \ml 



(28) 



The Boltzmann approximation can also be used to calcu- 
late the electronic contribution to other thermodynamic 
quantities in terms of the band gap and the effective 
masses. 



u{T = 0) + 



/ * \ 3/4 

[m^mi^) exp 



X (fcsT)^/' (SfefiT + Eg) , (29) 
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FIG. 10: (a) Effective masses of electrons in the first band gap. (b) Effective masses of holes in the first band gap. (c) The 
effective masses of electron in the second band gap. The dashed lines indicate the transitions from direct to indirect band gaps. 
The indirect band gaps are in between the two dashed lines, (d) The effective masses of holes in the second band gap. 
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(31) 



f-u{T = 0) 



X exp 



-Eg ' 

2kBT^ 



(32) 
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Here Eg is the band gap. The Bohzmann approximation 
assumes that there are no other kinks in the density of 
states besides the square root behavior described by Eq. 
(|27p within fc^T of the band edge. This approximation 
works best for small band gaps. For large band gaps the 
bands are narrow and there are Van Hove singularities 
near the band edges. The electronic contributions to the 
thermodynamic quantities of semiconductors and insula- 
tors are exponentially suppressed at low temperatures by 
the factor exp{—Eg/2kBT) and are often simply ignored. 
Equations can be used to estimate the tem- 

peratures where it is no longer reasonable to ignore the 
electronic contribution to the thermodynamic properties 
of semiconductors. 



VI. CONCLUSIONS 

The classic problem of electrons moving in a 1-D 
square-wave potential was considered and expressions 
were derived for the density of states and the effective 
masses of electrons and holes. These expressions allowed 
us to efficiently calculate the band structure and elec- 
tronic contribution to the thermodynamic properties in 



3-D separable square-wave potentials. This relatively 
simple model produces a wide range of band structures 
including metals, semiconductors, and insulators. Plots 
were presented showing the parameters for which band 
gaps appear. It was observed that for this class of po- 
tentials, band gaps only appear when V > 0.9Tr^h^ /ma^. 
Thus the condition for the existence of band gaps de- 
pends strongly on the lattice constant. The density of 
states at the Fermi energy and its derivative as well as the 
effective masses of the electrons and holes were also cal- 
culated. This makes it possible to estimate the thermo- 
dynamic properties using the standard approximations 
of the Sommerfeld expansion for metals and the Boltz- 
mann approximation for semiconductors. Expressions 
were derived for the thermodynamic properties of semi- 
conductors in the Boltzmann approximation. These sim- 
ple models assume that the density of states is a smooth 
function in an energy range fc^T wide. For those cases 
where a Van Hove singularity causes the density of states 
not to be a smooth function, the thermodynamic proper- 
ties were calculated numerically and the results are avail- 
able in the supplementary material. 
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